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Abstract 

The method analytic continuation of operators acting integer n-times to complex 
s-times (hep-th/9707206) is applied to an operator that generates Bernoulli numbers 
B n (Math. Mag., 70(1), 51 (1997)). B n and Bernoulli polynomials B n (s) are 
analytic continued to B(s) and B s {z). A new formula for the Riemann zeta function 
C(s) in terms of nested series of C( n ) 1S derived. The new concept of dynamics 
of the zeros of analytic continued polynomials is introduced, and an interesting 
phenonmenon of 'scatterings' of the zeros of B s (z) is observed. 



1 Introduction: Bernoulli Numbers 



Bernoulli numbers B n were discovered by Jakob Bernoulli (1654-1705) Q. They 
are denned 01 H as 



% = £-T* n > W<2tt, n = l,2,3,...GZ+ (1) 
- 1 * — ' n! 

n=0 



Expanding the l.h.s. as a series and matching the coefficients on both sides gives 

Si = -1/2, B n {zl ' ° ddn (2) 
1 ' "1^0, even n v ' 

With this result, (0) can be rewritten as 



oo „ 



e 2 - 1 2 ^ 2n 

71= 



Alternatively, Bernoulli numbers can be defined as satisfying the recurrence relation 



n—l 



Bernoulli numbers are interesting numbers. They appear in connection with 
a wide variety of areas, from Euler-Maclaurin Summation formula in Analysis || 

and the Riemann zeta function in Number Theory || |?|, to Rummer's regular 
primes in special cases of Fermat's Last Theorem and Combinatorics H. 



2 A Tree for Generating Bernoulli Numbers 

It was shown in [Q how a binary Tree for generating Bernoulli numbers can be 
intuited step-by-step and eventually discovered. In the process of calculating the 
analytic continuation of the Riemann zeta function to the negative half plane term- 
by-term, an emerging pattern was observed. The big picture of the structure of 
the Tree became apparent on comparing the derived expressions with the Euler- 
Maclaurin Summation formula. 

In this paper, we start with the Tree and proceed on to find interesting appli- 
cations. While doing so, we will encounter some surprising consequences. 

The Tree can be constructed using two operators, Ol and Or. 

At each node of the Tree sits a formal expression of the form — — - . 

a\b\ . . . 

Define Ol and Or to act only on formal expressions of this form at the nodes 
of the Tree as follows: 

° L : ^h.^ ( a + T)!&!... (5) 
±1 ±1 

° R ■ ^ja-. (6) 

Schematically, 

• Ol acting on a node of the Tree generates a branch downwards to the left 
(hence the subscript l in Ol) with a new node at the end of the branch. 
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Figure 1: The binary Tree that generates Bernoulli numbers. 



Or acting on the same node generates a branch downwards to the right. 



Form a finite series out of the sum of the two non-commuting operators 

/ n-l N 



S n = (0 L + OrT 



+ 1 

2! 



o n L + Y J OT 1 - k o R o k L + --- + o n R 



fc=0 



2! 



(7) 



This is equivalent to summing terms on the n-th row of nodes across the Tree. 
Bernoulli numbers are then simply given by 

B n = n! S n -i for n > 2 (8) 



eg., B 3 =3\S 2 = 3! (O l + Or) 2 



= 3! (0 L + Or) {O l + Or) 
+1 



+ 1 
2! 



3! (O l O l + L R + OrO l + OrOr) ( — 



= 3! 



+1 ^ , ^ +1 \ 
4! + 2!3! + 3!2! + 2!2!2! 







By observation, this Sum-across-the-Tree representation of S n is exactly equiv- 
alent to the following determinant known to generate B ni 
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(9) 
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3 Analytic Continuations 



3.1 Analytic Continuation of Operator 



First, we introduce the idea of analytic continuing the action of an operator following 
0]. We are used to thinking of an operator acting once, twice, three times, and 
so on. However, an operator acting integer times can be analytic continued to an 
operator acting complex times by making the following observation: 

A generic operator A acting complex s-times can be formally expanded into a 
series as 

1 



A s = Cwl-[wl-A]^ =w s (l 



\ n— 1 

t oo 

v n— 1 



-i) n 



-i) n 



w 

n-1 

n(*-*) 

k=Q 
n—1 

H(s-k) 

lk=0 



1 



-.4 



-.4 



A' : 



(10) 



where s G C, w € R, and 1 is the identity operator. 



The region of convergence in s and the rate of convergence of the series will in 
general be dependent on operator A, parameter w, and the operand on which A 
acts. 



3.2 Analytic Continuation of the Tree-Generating Operator 



Just as in (|10[), the Tree-generating operator (Ol + Or) acting (s — 1) times on 
-J— - J can be similarly expanded as 



2! 



-1 



/ oo 



w 



(-1)" 



l[(s-k) 



U-l 



l 7/1 



1 , y(~l) T 

2 n\ 

\ n—1 



Y[( S -k) 



Lfc=l 



1 A /-l 

9 + E — 



rn—1 



B m _|_i 
mj (m+1)! 



which converges for Re(s) > (1/w) where s G C, i»eR, w > 0. 



(ii) 



3.3 Analytic Continuation of Bernoulli Numbers 



-1 



B n = nl (O l + On)-' ( — ) = r(l + n) (0 L + Or)*- 1 ( — 



-1 



2! 



(12) 



Now that we can analytic continue the Tree-generating operator (Ol + Or) with 
(11), if we do so, we turn the sequence of Bernoulli numbers B n into their analytic 
continuation — a function 



B(s) = T(l + a) {Ol+O r ) 



■- M ±1 
2! 



(13) 



(~l) r 



U(s-k) 



,fc=i 



i M 



n\ B 



m+l 



my (m+l)! 
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which converges for Re(s) > (1/w) , real w > 0. 

So effectively, by the method of analytic continuation of operator, we have now 
obtained the function B(s) as the analytic continuation of Bernoulli numbers. 



B(s) 




lO 12 14 



Figure 2: The curve B(s) runs through the points of all B n except B\. 



All the Bernoulli numbers B n agree with -B(n), the analytic continuation of 
Bernoulli numbers evaluated at n. 



except 



B(n) = B n for n > 2 
B(l) 



- but Bi = - \ 
2 2 



(14) 
(15) 



4 Missing Signs in the Definition of B n 

Looking back at (jl|) to we can see that the sign convention of B\ was actually 
arbitrary. ( fLq ) suggests that consistent definition of Bernoulli numbers should really 
have been 



B„ 



t = Ec- 1 )"^ 

1 ^— ' n! 



|z| < 2vr , n = 1,2,3, ... e 



n=0 



fe=0 



n + 1 



Bit 



Sn = 1 



(16) 



(17) 



which only changes the sign in the conventional definition of the only non-zero odd 
Bernoulli numbers, B x , from B x = -1/2 to B x = B(l) = 1/2 . 

So here's my little appeal to the Mathematics, Physics, Engineering, and Com- 
puting communities to introduce the missing signs into the sum in the definition of 
Bernoulli numbers as in (16) and (^j) because the analytic continuation of Bernoulli 
numbers fixes the arbitrariness of the sign convention of B\. 



5 A New Formula for the Riemann Zeta Function 



Bernoulli numbers are related to the Riemann zeta function as 

B n +i 



C(-n) 



n + 1 



(18) 



4 



sv ; 2(2n)! 

for n = 0,l,2,...e Z+. 

From the above analytic continuation of Bernoulli numbers, 

B n i * 



(19) 



C(-n) = - 



71+1 



=> C(l-s) = - 



B{s + 1) 

8+1 



(20) 



Replacing B(s) in ( J2C| ) with the series in (13) and noting that T(l + s)/s = T(s) 
gives 



C(l-s) = - u/ 



(-1)" 



Lfe=l 



-1Y / n\ B 



>m+l 



— ' V w J \mj (m +1)! 

(21) 

which converges for Re(l — s) < 1 — {1/w) , real w > 0. 

The functional equation of the Riemann zeta function relates £(1 — s) to £(s) as 



C(i- S ) = 2(27r)- S r( s ) COs(^) CW 

Applying this relation to ( pTl) yields 



(22) 



7T5 

cos ( y ) C(s) 



9 9 2-^ «! 



n=l 



.fc=i 



l\/n\ B 



>m+l 



V w ) \ m J ( m +1)! 

(23) 



or in the limiting form 



C(a) = — hm { 



(-1)" 



2 + ^ ;/! 

\ n— 1 



H(S-k) 



Lk=l 



— 1 \ / n \ B, 



^.\w J \m J (m + 1 ) ! 



>m+l 



COS 



(-1)" 



_^ OO 

2 + ^ ,/! 

n— 1 



IKS-*) 



5-E 



y ml 



2 * — ' \ w / \mj m 

m— 1 



(2t™) s 
— hm < 



(24) 

a nested sum of the Riemann zeta function itself evaluated at negative integers, 
which converges for Re(s) > (1/w) , real w > , and the limit only needs to be 
taken when s = 1, 3, 5, . . . S Z^ dd , the set of positive odd integers, for which the 

/7TS\ 

denominator cos J = 0. 



This is consistent with 
(-l)" +1 (27r) 2 ™ 



C(2n 



2 (2n)\ 



Bn 



C(2s) 



{ 2«r r B( 2.) 



2T(l + 2s) Lcos(tts) 



COO 



lim 



2T(l + s) s^S | /7TS 

cos ( - 



(25) 
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Figure 3: W(s) 




envelopes ± \C(s)\ shown as dashed curves. 



Since both B(k) and cos ( — j = for odd integer k € Z^ dd \{l}, L'Hopital 



rule can be applied to the limit giving 



lim 

s—>s 



B(s) 



B(s) 



(t) 



= < 



s i z; 



s = 1 



(26) 



n/2 2B'(«) 



. »ez+„\{i} 



where the prime denotes differentiation and so B'(s) is the derivative of the function 

£(*). 

Therefore we have now found the apparently missing 'odd' expression 'dual' to 
the 'even' expression (|l9|). 



(even) : 



(odd): 



-i(2tt) 2 " B 2 



C(2n) = (-l)" +i ^TTT 
C(2n+1) = (-1) 



(2n)! 2 
„(2^) 2 " +1 B'(2n + 1) 



(271+1)! 7T 

for n = 1,2,3, ... G Z+. 

In fact, we can express -B'(s) in terms of C'( s )> the derivative of £(s). 

= -sC(l-s) 
B'(s) = -C(l-s) + sC'(l-s) 
S'(2n+1) = (2n+ 1) C'(-2n) 

as C(-2n) = 

1)n (2^) 2 " +1 



(27) 
(28) 



(29) 



C(2n + 1) = (-1)' 



(2n)\ 7T 



C'(-2n) for n e Z+. (30) 



This is just the differential of the functional equation and so verifies the consis- 
tency of B(s) and B'(s) with B n and ((s). 
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6 Other Undiscovered Half of Bernoulli Numbers 

From the relation (p0|), we can define the other analytic continued half of Bernoulli 
Numbers 

B(s) = -sC(l-s), B(-s) = s((s + l) 
=S> B_ n = B{-n) = n C(n + 1) , n G Z+ (31) 
Since £(n + 1) — > 1 as n — > oo, ~ n asymptotically for (—ri) <C —1. 




Figure 4: The curve -B(s) runs through the points and grows ~ n asymptot- 
ically as (—n) —* — oo. 



7 More Related Analytic Continuations 

7.1 Analytic Continuation of Bernoulli Polynomials 

The conventional definition (2| j| of Bernoulli polynomials B n (x) also has an ar- 
bitrariness in the sign convention. For consistency with the redefinition of B n — 
B(n) in ( |l6| ) and (0), Bernoulli polynomials should be analogously redefined as 



(iz)' 
. . , . ml 

\k=0 / \m=0 / n=0 



=* B„(x) = 1 B k x-- k (32) 

fc=0 ^ ' 

The analytic continuation can be then obtained as 

n h sGR, s>l, a: i— > u; £ C 

B k i ^ + («-[*])) = -(fc +(*-[*])) C(l-(fc + (s-H))) 
'n\ T(l + s) 



v fcy r(i + * + («-[«])) r(i + [a] -fe) 

W (l)n +fc r(1+g) B (k+{ S -[S])) 

=> Bn{x) » b( S , w) = j: T{l+k+is _ W)) r(1 + [s] _ k) 
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M+i (1) n +fc+ i r(1 + fi) B ( (fc _ 1) + (s _ [a]) ) w [s ]+ i-k 

£f r(* + («-[«])) r(2 + [*]-fe) 1 j 

where [s] gives the integer part of s , and so (s— [s]) gives the fractional part, 

eg., B 2 (x) = 5(2, x) = 1/6 - x + x 2 w 0.16667 - x + x 2 

B(2.01,x) w 0.16420 - 0.99660x4- 1.00920 x 2 - 0.00554 a; 3 
B(2.99,x) w 0.00092 - 0.50576 x + 1.50001 x 2 - 0.98744 x 3 
B 3 (x) = B(3,x) = -0.5 x + 1.5 x 2 - a; 3 



B(s,x) o 



^^^^^^^^ 


^ 2.2J,.\-; 







— 1 O 1 2 



Figure 5: Deformation of the curve of B?,{x) — B(2,x) into the curve of B^(x) — 
B(3,x) via the real analytic continuation B(s,x), 2 < s < 3, i£ R. 



7.2 Analytic Continuation of Euler Numbers and Euler Poly- 
nomials 



Euler numbers E n and Euler polynomials E n (x) are defined |2J [|3j as 

2 e z z n 

sech z = — = } E n — , \z\ < 2tt , n e Z + ; # n = for odd n (34) 



71=0 



2 e xz 2 e z / 2 e 2 ^" 1 / 2 ) / v-^ 



CO I. 

Z 



„2z+l 



2 k k\ 1 



~ (x - \) m z m 

< 4 rn I 



m=0 



n=0 



fc=0 



and are related to Bernoulli polynomials as 

2 



E n (x) = — (B n+1 (x) - 2" +1 B n+1 (f ) 



2 n+1 B, 



n+l 



z — ■ n! 

(35) 

(36) 
(37) 
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Their analytic continuation then follows straightforwardly from the analytic con- 
tinuation of Bernoulli polynomials 

B n {x) h-> B(s,w) , E n >-^E(s), E n {x)^E{s,w) 

where s £ R , s > 1 . ihkieC. 

8 Beautiful zeros of Bernoulli Polynomials 

8.1 Distribution and Structure of the zeros 

Zeros of Bernoulli polynomials are solutions of B n (w) — 0, w G C. 



Im(w) 



-7.5 -5 -2.5 2.5 5 7.5 10 

Re(w) 



Figure 6: Zeros of Bernoulli polynomials Bso(w). 
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Figure 7: Stacks of zeros of Bernoulli polynomials B n {w) form a 3-D structure. 



The real zeros, except the outermost pair in general, are almost regularly spaced, 
while the complex zeros lie on arcs symmetrical about Re(w) = 1/2. 
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8.2 Observations, Theorem, Conjectures, and 
Open Problems 

1. Symmetries 

Prove that B n (w), w G C, has Re(iu) = 1/2 reflection symmetry in addition 
to the usual \m(w) — reflection symmetry analytic complex functions. The 
obvious corollary is that the zeros of B n (w) will also inherit these symmetries. 

If B n (w Q ) = 0, then B n (l - w ) = = B n (w f ) = B n (l - tuj) (38) 
where ' denotes complex conjugation. 



2. N on- degenerate zeros 

Prove that B n (w) = has n distinct solutions, ie., all the zeros are non- 
degenerate. 



3. Central zero Theorem 

If-B2 ro +i nas R e ( w ) = 1/2 an d \m(w) — reflection symmetries, and (2m+l) 
non-degenerate zeros, then (2m) of the distinct zeros will satisfy (p8|). If the 
remaining one zero is to satisfy (|38| ) too, it must reflect into itself, and therefore 
it must lie at 1/2, the center of the structure of the zeros, ie., 

B n (l/2) = V odd n (39) 

4. Counting of real and complex zeros 

Prove that the number of complex zeros Zc(n) of B n (w) lying on the 4 sets 
of arcs off the real plane, \m(w) ^ 0, is 



z C(r 



, where 5 ■■ 



n + 30 
21 



(40) 



and [ ] denotes taking the integer part. The factor 4 comes from the above 
2 reflection symmetries. 

Since n is the degree of the polynomial B n (w), the number of real zeros ZR( n ) 
lying on the real plane \m(w) = is then z R( -„) = (n— zc(n))- 

See Appendix for tabulated values of zr(„) and Zc(n)- 

5. Asymptotic Lattice of real zeros 
It is known that from (|7.l|), 



J2 {B n (w + 1) - B n (w)) - = ze wz = ^ 7 Try 

n=0 n=l v ' 

=> B„(«i + 1)-B„H = u;™- 1 

S„(l) = S n (0) = B n = V odd n > 3 

Show that all the real zeros of B n (w) except the outermost pair in general, 
are approximately regularly spaced at the staggered lattice points 

w _ f 0, ±1/2, ±1, ±3/2,... eZ/2 for odd n ^ 

\ ±1/4, ±3/4, ±5/4, ... eZ/2 + 1/4 for even n 
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Re(w) 



Figure 8: Inner real zeros converge to a staggered Lattice structure. 



and becomes increasingly located exactly at these lattice points as n — > oo. 
See Appendix for tabulated solutions of B n (w) = . 

6. Relation between zeros of Bernoulli and Euler polynomials 

Choose any zero of Bernoulli polynomial B n {w) and denote it as Ws( n ), 
ie. B n (w B(n) ) = 0. 

Prove that 

asn^M, E n (2 (w B[n) - 1/2) + 1/2) = E n (2 w B{n) - 1/2) = (42) 

ic., the structure of the zeros of Euler polynomials resembles the structure of 
the zeros of Bernoulli polynomials but doubled in size in the limit the degree 
of the polynomials n — > oo. Both structures are centered at w = 1/2. 

7. Bounding Envelopes and Trajectories of complex zeros 

Find the equation of envelope curves bounding the real zeros lying on the 
plane, and the equation of a trajectory curve running through the complex 
zeros on any one of the arcs. 

8.3 Dynamics of the zeros from Analytic Continuation 

Bernoulli polynomial B n (w) is a polynomial of degree n. Thus, B n (w) has n zeros 
and B n+ i(w) has (n+ 1) zeros. When discrete n is analytic continued to continuous 
parameter s, it naturally leads to the question: 

How does -B(s,w), the analytic continuation of B n (w), pick up an additional 
zero as s increases continuously by one? 

This introduces the exciting concept of the dynamics of the zeros of analytic 
continued polynomials — the idea of looking at how the zeros move about in the 
w complex plane as we vary the parameter s . 

Continuity shows that the additional zero simply "flows in from infinity" . 
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To have a physical picture of the motion of the zeros in the complex w plane, 
imagine that each time as s increases gradually and continuously by one, an addi- 
tional real zero flies in from positive infinity along the real positive axis, gradually 
slowing down as if "it is flying through a viscous medium" . 

For s < 5, the additional zero simply joins onto the almost regularly spaced 
array of the real zeros streaming slowly towards the negative real direction. The 
array of zeros continue to drift freely until one by one they hit a bounding envelope 
which grows in size with s. 

As s approaches every integer > 5, an interesting phenomenon occurs: A pair 
of real zeros may meet and become a doubly degenerate real zero at a point, and 
then bifurcate into a pair of complex zeros conjugate to each other. Thus, the pair 
of real zeros appears to "collide head-on and scatter perpendicularly" into a pair of 
complex zeros. 




Re(w) 



Figure 9: The flow of zeros of B(s, w) forms a complex 3-D structure. 
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3 fundamental kinds of scattering can be observed: 




(a) Point (b) Loop (c) Long-range 



Figure 10: 3 fundamental kinds of scattering. 



• Point scattering: 

A pair of real zeros scatter at a point into a pair of complex zeros which head 
away from each other indefinitely. 

• Loop scattering: 

The same as point scattering but the pair of complex zeros loops back to 
recombine into degenerate real zeros within unit interval in s and then scatter 
back into a pair of real zeros, much like the picture of pair production and 
annihilation of virtual particles. 

• Long-range sideways scattering: 

The additional zero that flies in appears as if "it is carrying with it a line 
front of Shockwave" that stretches parallelly to the \m(w) axis. When the 
"shockwave" meets a pair of complex zeros that are looping back, the pair 
gets deflected away from each other momentarily before looping back again, 
while the additional zero gets perturbed and slows down discontinuously. 

The whole complex structure can then be reduced to simple combinations of 
these 3 kinds of scattering. 




Figure 11: A composite of Point, Loop and Long-range scatterings. 
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The movies (animated gifs) showing the motion of zeros in the complex w plane 
at slices of s can be downloaded from Internet at |ll| and viewed with Netscape or 
Internet Explorer WWW browser. 

8.4 Open Challenges in Generalization 

1. Generalize the above results of B(s,w) consistently from s £ R , s > 1 to 
s G C. 

2. Derive a set of expressions which give the values of s where the point, loop 
and long-range scatterings occur. 

9 Conclusion 

Bernoulli numbers and polynomials appear in many areas. In particular, if we 
assume the Riemann Hypothesis, when s G C, ( pif ) should converge to zero only on 
the line s — 1/2 ± it. This remains to be proved. More of this aspect is analysed 
in©. 

To further self-explore these fascinating properties, feel free to download and 
adapt the executable Mathematica codes from fl2|| . Have fun! 

In the meantime, it would be interesting to imagine what Jakob Bernoulli and 
Euler would say on these analytic continuations of their numbers and polynomials. 



Acknowledgement 

Special thanks to Y.L. Loh, W. Ballman, B. Lui, P. D'Eath, and K. Odagiri for 
discussion, all the friends in Cambridge for encouragement, and Trinity College and 
UK Committee of Vice- Chancellors and Principals (CVCP) for financial support. 



References 

J. Bernoulli, "Ars Conjectandi" 

H. Bateman, Higher Transcendental Functions, Vol 1., (McGraw-Hill, 1953). 

M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, (Dover, 
1970). 

S.C. Woon, Math. Mag. 70(1), 51 (1997). 

M. Spivak, Calculus, (Benjamin, 1967), 482, Problem 17. 

E.C. Titchmarsh, The Theory of the Riemann zeta-function, (Oxford, 1986). 

S.C. Woon, Chaos Solitons & Fractals 5(1), 125 (1995). 

P. Ribenboim, The little book of Big Primes, (Springer- Verlag, 1991). 

S.C. Woon, "An alytic Continuati on of Operators — Operators acting complex 
s-times", e-Print licp-th/9707206 . 



[10] S.C. Woon, "Chaos, Order, and 2 Constants in the Riemann Zeta Function" , 
(e-Print archive chao-dyn, in preparation). 



14 



[11] S.C. Woon, Movies of scattering zeros, (1997) 

http: / /www. damtp.cam.ac.uk/user/scw21/papers/97051/movie. html 

[12] S.C. Woon, Mathematica codes, (1997) 

http:/ /www. damtp.cam.ac.uk/uscr/scw21/papers/97051/codes. html 



Appendix: 

Number of real and complex zeros of B n (w) 



B n (x) 
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complex 
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16 


26 
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B n {x) 
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complex 


degree 
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ii 


z R(n) 


z C(n) 
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24 
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41 
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43 


11 


32 


44 


12 


32 


45 


13 


32 


46 


14 


32 


47 


15 


32 


48 


12 


36 


49 


13 


36 


50 


14 


36 


51 


15 


36 


52 


16 


36 


53 


17 


36 



B n {x) 


real 


complex 


degree 


zeros 


zeros 


n 


ZR(n) 


zc(n) 


54 


14 


40 


55 


15 


40 


56 


16 


40 


57 


17 


40 


58 


18 


40 


59 


15 


44 


60 


16 


44 


61 


17 


44 


62 


18 


44 


63 


19 


44 


64 


16 


48 


65 


17 


48 


66 


18 


48 


67 


19 


48 


68 


20 


48 


69 


17 


52 


70 


18 


52 


71 


19 


52 


72 


20 


52 


73 


21 


52 


74 


22 


52 


75 


19 


56 


76 


20 


56 


77 


21 


56 


78 


22 


56 


79 


23 


56 


80 


20 


60 



15 



Solutions of B n (w) = 



odd n 



w 



7/3 



-0.26376, - + 



2) 



1 V9T2V37 
2\/3 



1.26376 



1 V9 - 2V3? 

2 2V3 

1 ^9 - 2V3« 

2 + 2^5 



1 y / 9 + 2y / 3~ i 

2 + 2^ 



-0.38137- 0.16376 i, 
-0.38137 + 0.16376 i, 
1.38137- 0.16376 i, 
1.38137 + 0.16376 i 



even n 



w 



2V3 



V15 


f 2 V 


/30 


2i 


/15 




\/l5 


-2 V 


/30 


2i 


/15 




Vl5 


-2 V 


/30 


2i 


/15 




\/l5 


f 2v 


/30 


2i 


/15 





0.21133, ( J + ^7= 
1 2 2V3 



-0.15770, 
0.24034, 
0.75967, 
1.15770 



0.78868 



16 



